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Two-body matrix elements of arbitrary local interactions are written as the sum of separable terms 
in a way that is well suited for the exchange and pairing channels present in mean-field calculations. 
The expansion relies on the transformation to center of mass and relative coordinate (in the spirit 
of Talmi's method) and therefore it is only useful (finite number of expansion terms) for harmonic 
oscillator single particle states. The converge of the expansion with the number of terms retained is 
studied for a Gaussian two body interaction. The limit of a contact (delta) force is also considered. 
Ways to handle the general case are also discussed. 
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I. INTRODUCTION 

The evaluation of the pairing field in mean field the- 
ories like Hartree- Fock plus Bardeen-Cooper-Schrieffer 
(HF+BCS) or Hartree-Fock- Bogoliubov (HFB) is a com- 
putationally intensive task due to the non-local character 
of the pairing tensor and the effort is comparable to the 
one devoted to the evaluation of the exchange field in 
the Hartree Fock method. Therefore, zero range pairing 
interactions are thoroughly used in standard mean field 
calculations in order to reduce the computational bur- 
den (see Ref [l| for a recent review). The price to pay 
for the use of zero range pairing forces is the introduc- 
tion of an "active space" around the Fermi level to cut 
away the ultraviolet divergences inherent to any contact 
interaction (the pairing matrix elements are independent 
of the momentum transfer in nuclear matter). Thus, it is 
customary to take into account in the gap equation only 
those single particle levels lying inside a so called "active 
window'' around the Fermi level and whose (sharp or soft) 
boundaries are defined using reasonable (but arbitrary) 
assumptions. The boundaries of the window as well as 
the pairing interaction strength are usually fixed locally 
but they are usually not allowed to vary as a function 
of other relevant degrees of freedom like the quadrupole 
deformation of the nucleus. As a consequence, the im- 
pact of the rigid definition of the "pairing active window" 
in some observable magnitudes can be relevant (see, for 
instance, [H for a discussion on the impact of the win- 
dow's size on fission barriers). From this perspective, the 
use of finite range pairing forces seems to be quite un- 
avoidable and this is a common argument to praise the 
use of Gogny type interactions Q in HFB or HF+BCS 
like mean field calculations |4j . It has to be mentioned 
that a renormalization scheme for the phenomenological 
contact pairing interactions has been proposed in Refs 
0, [|| to cure the "active window" problem. Recently, 
proposals to use as pairing interaction in finite nuclei a 
realistic finite range two body bare interaction Q or a 
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low momentum evolved version of it [7-9] have been dis- 
cussed in the literature: the essence of the proposals is to 
introduce a separable interaction in momentum represen- 
tation [7J that is used to fit the realistic interaction (or 
a low momentum evolved version of it) by resorting to 
fitting protocols involving nucleon-nucleon phase shifts 
or pairing gaps in nuclear matter. As a consequence of 
the separable character in momentum space the force is 
non-local in real space but in a way that simplifies the 
evaluation of the matrix elements needed for finite nuclei 
calculations [1, @. Separable interactions were also used 
as a way to tackle the solution of the Lippman-Schwinger 
equation for the T operator with realistic nucleon-nucleon 
potentials [10]. Recently, a separable form in momentum 
space of the pairing interaction of a Gaussian two body 
force (as in the Gogny force) has been proposed [IlT - fl3j 
as an alternative to more standard approaches based on 
zero range contact pairing interactions. As the approxi- 
mation relies on the same ideas of Refs. 0,13 a non-local 
form of the approximate interaction in real space is ob- 
tained. The non locality of the interaction suggests the 
introduction of the Talmi-Brody-Moshinsky transforma- 
tion [T^ - tl^ ] to center of mass and relative coordinate to 
simplify the evaluation of the pairing matrix elements 
in the harmonic oscillator basis. Explicit expression for 
the matrix elements are given in |l ll | for an harmonic 
oscillator basis with spherical quantum numbers and a 
Gaussian interaction. Next, it has been shown that this 
new form is consistent with the RPA framework [12| . Fi- 
nally, the scheme has been extended [13] to the case of 
a harmonic oscillator basis with axial quantum numbers. 
The extension to the case of a harmonic oscillator basis 
with triaxial quantum numbers is straightforward. In the 
present article I show that the transformation to center 
of mass and relative coordinates can be used for general 
two body interactions to obtain a kind of separable ex- 
pression for the pairing two body matrix elements which 
involves the one body matrix elements of the interac- 
tion for the relative coordinate wave functions. The use 
of the spectral representation of this interaction matrix 
leads to a more exp licit separable form of the kind con- 
sidered in 0-0 Illl4l3l | . By using two specific examples 
in one dimension it is shown that the number of sepa- 
rable terms in the expansion can be severely cut down 
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for short range interactions. The validity of the expan- 
sion concerning ultraviolet divergences is also discussed. 
Finally, computational schemes to deal with the Yukawa 
(and Coulomb!) potential are discussed. It has to be 
stressed that the present approach is not limited to pair- 
ing matrix elements and can easily be extended to deal 
with the exchange matrix elements required for the evalu- 
ation of the Fock potential and therefore it could be used 
as an approximation scheme to shorten the numerical 
burden of HFB calculations with finite range forces like 
Gogny 0||, other recent proposals based on the Yukawa 
potential [T3| and even the Coulomb potential. 



© is relevant. Explicit forms of the TBMC coefficients 
have been discussed several times in the literature both 
in the 3D spherical form of the HO wave functions [14- 
[H| as in the one-dimensional case [2(| . Introducing the 
expansion of Eq. ^ into the definition of the matrix 
element we obtain the result 



x ni,ri2,n 3 ,n4 



= £E M »£>»' M ££' ( 2 ) 



/V 



where we have made use of the orthogonality of the center 
of mass wave functions <j>pf(R) and introduced the matrix 
elements 



II. SEPARABLE APPROXIMATION TO 
TWO-BODY MATRIX ELEMENTS 

A. General procedure 

Let us consider the (not antisymmetrized) 
two-body matrix element ^n 1 .n2,n 3 .n 4 = 
Jdri/df20* 1 (f ! i)^ 3 (f2)i;(fi - r 2 )0„ 3 (ri)(/>„ 4 (r^) 
for harmonic oscillator (HO) wave functions (j) n (r) 
(this is not a fundamental restriction as more general 
single particle wave functions can always be expanded 
in a HO basis). Taking into account that at the end 
we will be folding this matrix elements with a matrix 
(either the pairing tensor or the density for the exchange 
potential) with indexes 713 and 714 the best approach 
for the evaluation of this matrix element is that of 
Talmi [TJ| and Brody-Moshinsky [H [TJ| (see also 
[H . [HI for a general discussion). Within this method, 
the product of HO single particle wave functions with 
different arguments 0„ 3 (r\)<f> ni (r^) is expanded in terms 
of suitable center of mass R = r^(^i + ft) an d relative 
coordinate r= ^(^i — ft) wave functions (this unusual 
definition of the center of mass and relative coordinate 
renders some of the expressions simpler than with the 
usual definition) 

K 3 in)<t> ni (ft) = MZl^ N (R)Mr1 (1) 

Nn 

The center of mass and relative coordinate wave func- 
tions have the same structural form as the original wave 
functions and the expansion is finite (the range of val- 
ues of N and n is finite !). This is a direct consequence 
of the special structure of the HO wave functions (the 
product of a Gaussian times a polynomial). Apart from 
the HO, this peculiarity is only preserved in the case of 
plane waves and the developments considered below can 
be extended easily to this case too. The expansion coef- 
ficients M^™ 2 are referred to as Talmi-Brody-Moshinsky 
[TM3| coefficients (TBMC). The TBMC have a selection 
rule that will help to reduce the number of final separable 
terms, it is n\ + n 2 = N + n (see Appendix A for the gen- 
eral expression and selection rules in the one-dimensional 
case) which implies that only one of the two sums in Eq. 



Vnn' = / df4>l(r)v(V2r)(j) n '(r). 



(3) 



Please note the V2 factor in the argument of the interac- 
tion that is a direct consequence of the definition of the 
relative coordinate. This result can be easily generalized 
to non-local interactions of the form 



(rir 2 |u|ri^) = S(R - R')vftr*) 



(4) 



where R and r are the center of mass and relative coor- 
dinates associated to f[ and r^- Applying the transfor- 
mation of Eq. (fTJ) we obtain for the HO matrix elements 
the same expression as in Eq. © but replacing v nn i of 
Eq. © by 

»W= I ' dr I ' dr'<t>l(r)v{r,r')cl> nl {r'). (5) 



The v nn i are matrix elements of a hermitian matrix 
and therefore can be written, by resorting to the spectral 
decomposition of the matrix, as 



, = V n* 

J nn / 
L 



D* nL v L D n , L 



(G) 



where the coefficients vl (the eigenvalues of v nn > ) are real 
quantities. By introducing the coefficients 



(7) 



we can finally cast Eq. @| in a form that corresponds 
clearly to a separable expansion for pairing and exchange 
matrix elements as 



N L 



(8) 



This separable expansion is still exact and can be used to 
express the pairing field also as a sum of separable terms 

Aniri 2 = ^ l/ n 1 ,n 2 ,n 3 ,n4, K n 3 n4 — ^ ^n\n->. V L^NL 
ri3«4 NL 

with A NL = J2n 3ni ^n 3 n 4 K n 3 n 4 - At this point it is worth 
mentioning that this procedure can be straightforwardly 
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extended to the evaluation of the exchange field by in- 
troducing the new TBMC as 

K (fMn* ft) = E M^ N (R)Mrl (9) 

JVn 

to obtain 

— M Nn „ M* Nn ' 

ty n 1 ,n 2 ,n 3 ,n 4 — / J / J m n\m v nn ' m n 2 n 3 

N nn' 

which can be used for the evaluation of the exchange 
field. This possibility has not been considered in Refs. 

mm. 

The drawback of the (still exact) separable expansion 
being considered is that the number of terms involved is 
enormous unless some property of the interaction is in- 
voked to restrict it. To do so, we can use the general ideas 
of linear algebra on how to approximate an arbitrary ma- 
trix in terms of low rank Kronecker tensor products (see 
[2ll . |22| for an introduction to the subject) to justify the 
rank one approximation 

V ~ d *n d n> (10) 

where d n — ^/voD n o and we are assuming that the L = 
term in Eq. © corresponds to the largest eigenvalue 
Do of the interaction matrix v nn i. In this way we have 
reduced the number of separable terms by a factor that 
corresponds to the dimension of the matrix v nn i . For bet- 
ter accuracy, we can pursue further this idea to suggest 
the approximation 

L C 

Vnn' ~ 2^ D* nL V L D n > L (H) 

(where vl with L = 0, . . . , Lq stand for the Lq + 1 
largest eigenvalues of the matrix v nn >). Obviously, low 
rank Kronecker tensor products do not necessarily have 
to be based on the eigenvectors of the interaction ma- 
trix v and other alternatives inspired in the physics to 
be described could be incorporated easily (a possibility 
that comes up immediately is to use for d n in Eq. (|10[) a 
linear combination of the eigenvectors of v with weights 
physically inspired by the problem to be treated and to 
be fitted for the optimal reproduction of some quantity.) 

B. Zero range interaction 

In the following we will consider some typical interac- 
tions to asses the validity of the preceding approximation 
scheme and we will start by considering a zero range con- 
tact interaction v(r) — GS(r). In this case, the matrix 
v nn ' is simply given by 

iW = 257a&(0)MO) (12) 

which implies, in the language of Eq. ([5]), that all the 
eigenvalues vl are zero except the L = one. For this 



interaction, the rank one approximation of Eq. (|10p is 
exact and the reduction in the number of separable term 
is quite significant. Introducing the quantity 

n 

we finally obtain 

N 

Apart from the analytical result, this formula is telling 
us that for short range interactions the expansion of Eq. 
(Ilip is expected to be accurate enough by considering 
only a limited number of terms. 

In order to get a deeper insight into the number of sepa- 
rable terms in the preceding expansion for the zero range 
force we will focus on the one-dimensional case. Our ba- 
sis will be the one-dimensional harmonic oscillator basis 
4> n {x) with n — 0, . . . , M and containing M + 1 elements. 
The TBM expansion of the 1-D HO wave functions is an 
expansion of the product of two polynomials in the vari- 
ables x\ and X2 in terms of the product of polynomials in 
the variables X = +£2) and x = -h%{ x x ~ x "i)- Al- 

though its explicit expression is known since a long time 
[20l | . we present in Appendix A a brief derivation based 
on the generating function of the Hermite polynomials. 
As a consequence of the explicit form of the TBM coef- 
ficients, both quantum numbers N (center of mass) and 
n (relative) should range from up to 2M (i.e. 2M + 1 
possible values) if none of the selection rules of the TBM 
coefficients is taken into account. The parity selection 
rule (-l)«i+«2 = of the TBM coefficients can 

be used to reduce the number of terms in the N sum. As 
4> n (fl) is only different from zero for even values of n, the 
parity of TV has to be the one of n\ +ri2 and therefore the 
number of terms in the N sum of Eq. (I13[) gets reduced 
to half the initial value (M + 1 to be more precise). 

C. Gaussian interaction 

In order to analyze in more detail the consequences 
of the proposed expansion, we will study the exam- 
ple of a Gaussian interaction of range fi in one dimen- 
sion. The more realistic three-dimensional case can be 
elucidated from the results obtained here as, in that 
case, the matrix elements are the product of the one- 
dimensional ones along each of the three spatial dimen- 
sions. In addition, the zero range interaction results can 
be recovered in the limit y, — > providing further in- 
sight into the number of terms to be retained as well 
as potential problems with ultraviolet divergences. The 
one-dimensional HO wave functions will be denoted as 
(p n (z;b) = exp(-ifj-)(^„(z/6) where (p(z/b) is the poly- 
nomial part of the HO wave function given by the prod- 
uct of the Hermite polynomial of degree n times the nor- 
malization constant of the ID HO wave function. The 
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quantity to evaluate is 

f 2z 2 
v nn > = J dzipl(z;b)exp(-—^-)ip n >(z;b) 

= J dz<p*(z/b)exp(-(-^ + ^)z 2 )^ n ,(z/b). 

This integral can be carried out in many different ways 
but for our purposes it is better to use the transformation 
matrix D(b'/b) connecting ip n (z/b) with ip n (z/b') 

n 

(p n {z/b) = D nn ,(b'/b)£ n ,(z/b'), 
n'=0 

defined in Appendix B, to write 

v nn , =J2D* nL (B/b)D n , L (B/b) 

L 

with B/b = m/Vm 2 + 2fr 2 = V- It * s better to write the 
preceding result as 

v nn ,=J2K L (vW L+1 Dn> L (v) (14) 

L 

to make explicit the dependence in the range of the Gaus- 
sian /x. Although this expression is not an explicit power 
expansion in 77 due to the dependence of D on this pa- 
rameter, it is close to be so because because the ratio 
between the L + 2 term of the sum and the L one is given 
by ry 4 (l + ^ 2 /(26 2 )) 2 (n-L)(n'-L)/((L + 2)(L + l)) that 
shows a dominant 77 4 behavior. As r\ typically goes as fi/b 
we can conclude that Eq ([T4]) is a kind of expansion on 
the range of the interaction [i. For small ranges [i — > the 
parameter r\ tends to [ij (V2b) and therefore the most sig- 
nificant term in the sum is L = 0. The remaining terms 
decrease as /j 2L and justify to cut the expansion of Eq. 
(P3J) at some small Lc 

v%F = fl £> nL(v)v 2L+1 D n , L (r ] ). (15) 

L=0 

Using the connection between a Gaussian of vanishing 
range and the delta function 

5(x\ — X2) = lim — = 

it is possible to obtain Eq. fT2")) from Eq. (| 14[) by appro- 
priately taking the fj, — > limit. As can be shown easily, 
the [i — > limit of D n o(r)) (the only remaining term is 
L = 0) is given by 7r 1 / 4 6 1 / 2 i^„(0) as required. 

At this point it is worth noticing that neither the eigen- 
values nor the eigenvectors of the v nn ' matrix in the 
Gaussian case are given by r] 2L+1 or by the coefficients 
D n 'L(rj) above (see Eq. (fT"4]) ). In order to understand the 
differences, the largest eigenvalues obtained by numeri- 
cal diagonalization in the case N — 40 and for b = 2fm 
are plotted in Fig. [1] as a function of 77. We observe 



in the log-log plot that for small values of r\ the eigen- 
values have a linear behavior that corresponds to r/ 2L+1 
for L = 0, . . . recovering in this case the result of Eq. 
(H3J). However, for larger values of ij this is not the case 
indicating that the eigenvalues of v nn i depart from the 
quantities in Eq. (|14p . An interesting observation re- 
garding the behavior of the eigenvalues is that they are 
a decreasing function as fi decreases, justifying the idea 
that keeping only the largest eigenvalues is a kind of short 
range expansion of the matrix elements. Obviously, if 
all the terms are considered, both the expansion of Eq. 
(Q3J) and the general one of the spectral decomposition 
of the matrix v nn i give the exact answer. The differ- 
ences will show up in the approximate case when only 
a finite number of terms is retained. We have checked 
that the spectral decomposition always gives the best re- 
sults if quantified in terms of the mean square deviation 

a v = ^Jj2nn'( v nn' ~ v^?) 2 / (N + 1), but none of the 

approximate matrix elements v^F reproduce the exact 
ones, as it may happen when using Eq. (|15[) . As a conse- 
quence of the selection rules of the D coefficients involved, 
the exact values of v nn i for low values of n and n' can 
be obtained from Eq. (|15|) even for small values of Lq. 
This is not in contradiction with what it was said about 
ct y , as the approximation based on Eq. (|15[) provides, for 
large values of n and n', worse matrix elements than the 
ones obtained using the spectral decomposition of v n n' ■ 
The analysis below of the two body matrix elements has 
been carried out in parallel using both possibilities and 
it has been found that, although the spectral decompo- 
sition provides again a better overall approximation for 
the two body matrix elements, those with small quantum 
numbers (the most relevant for the physics of the system 
to be described) are slightly better reproduced by the ap- 
proximations based in Eq. (fl5|). For the above reasons 
we will only show results based on Eq. (|15l) . 

The result obtained in Eq. (fT5")) for short ranges \x 
suggests to approximate v nn i by keeping only the L = 
term. However, taking into account that the transforma- 
tion coefficients D{rf) preserve parity, we realize that the 
L = approximation will make automatically zero the 
v nn ' matrix elements with odd n. Therefore, we have to 
include both L = and L = 1 as the leading order. The 
next to leading order will correspond to take L = 2 and 
3 and so on. In order to test the convergence rate of such 
an approximation we have performed calculations com- 
paring the exact matrix elements of the one-dimensional 
Gaussian interaction with those obtained using Eq. © 
and restricting the sum in L from up to a parameter de- 
noted Lc as well as the sum in N from up to Nc- The 
size of the one-dimensional HO basis considered is of 21 
shells (n — 0, . . . , 20) and there are 97241 non zero ma- 
trix elements. The results for a Gaussian of range /i = 0.7 
fm are depicted in Fig. [2] On the left hand side panels 
all the matrix elements are plotted (y axis, approximate 
values, x axis exact ones) for different values of Nc and 
Lc- On the right hand side panels, the same kind of plot 
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-2 -1 

log 10 (T|) 

Figure 1: Largest eigenvalues Xi of the interaction matrix in 
the case of a one-dimensional Gaussian interaction of range 
/i. The eigenvalues shown are plotted as a function of the 
parameter rj defined in the text. Notice the log-log character 
of the plot. 



is depicted but this time only the matrix elements with 
ni, ri2, «3 and n 4 smaller than 7 are shown. The reasons 
are two: first that the realistic single particle orbitals are 
usually close to the HO wave functions with the same 
quantum number n and therefore and expansion on the 
HO basis usually requires the n — 2,n and n + 2 states for 
an reasonable representation. The second reason is that 
only a limited number of orbitals around the Fermi level 
usually play a role in the pairing properties and therefore 
only the quantum numbers n required to accommodate 
A particles are required. Usually seven major shells (the 
corresponding quantity in the three-dimensional case) are 
required to well accommodate of the order of 200 parti- 
cles. Therefore, in the left hand side panels there are 
many matrix elements that will contribute little to the 
pairing field and can be safely disregarded for the dis- 
cussion of the quality of the approximation. That is, 
the matrix elements depicted on the right hand side pan- 
els are supposed to be the most relevant for the pairing 
properties of the nucleus. By looking at the plot we ob- 
serve how the bigger the values of Nc and Lc are the 
better the approximation is (the points gather around 
the y = x line). We even observe how the iV c = 10 and 
Lc = 1 approximation, including eleven separable terms, 
is already quite reasonable for the reduced set of relevant 
matrix elements (right hand side panels). In Fig. [3] the 
same kind of plots as the ones of Fig. [5] are presented 
but this time for a Gaussian twice the previous range, 
i.e. n = 1.4 fm. We immediately realize the worsening of 
the approximation, what is consistent with the previous 
findings that this is a kind of "short range" expansion. In 
this Figure we observe that the convergence with N c is 
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Figure 2: Approximate matrix elements (y axis) for different 
values of Lc and Nc are plotted versus the exact values (x 
axis) for the one-dimensional Gaussian interaction with range 
H — 0.7fm and harmonic oscillator length b z =2.0 fm. On the 
left hand side panels, all matrix elements with m, ri2, ns and 
m ranging from up to 20 (i.e. a total of 21 4 matrix elements) 
are depicted. On the right hand side panels only those matrix 
elements with m, 712, ^3 and 714 smaller than 7 are depicted 
(see text for details). 



much faster than the one with Lc because already the 
Nc = 10 numbers compare quantitatively well with the 
ones of Nc — 40 whereas for Lc the results with Lc = 3 
are much better than the ones for Lc = 1. 

Finally, in Fig. [4] we present a more detailed study of 
the convergence with Nc in the case of [i — 1.4 fm keep- 
ing fixed Lc to the reasonable value of 3. We observe 
how decreasing the value of Nc from 10 to 8 to 6 de- 
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Figure 3: Same as Fig. [2] but for a range p, = 1.4fm of the 
Gaussian interaction. 

grades the quality of the approximation for the relevant 
matrix elements (right hand side panels) but this is not 
significant and it is quite likely that even Nq = 6 (twelve 
separable terms!) will provide already reasonable values 
of the pairing tensor in calculations with real nuclei. This 
point can only be tested in realistic calculations and work 
along this direction is in progress. 

It is also interesting to compare the approach of Refs. 
041) ITU - HI^ with the present one. In those references the 
two body Gaussian interaction is replaced by 

(nr 2 \v\^,fi) =GS(R- R')^2\ a P a (r)P a (r') (16) 

a 

where R and f are the center of mass and relative co- 
ordinate of fi and r*2- The functions P a (f) are taken 
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Figure 4: Comparison between exact (a; axis) and approxi- 
mate (y axis) pairing matrix elements of an one- dimensional 
Gaussian interaction with range fi — 1.4 fm as a function of 
Nc (see text for details as well as Fig. [2] for the meaning of 
the left hand side and right hand side panels). 



as Gaussian with widths adjusted to reproduce nucleon- 
nucleon phase shifts 0, Q or the nuclear matter pairing 
gap (nHl3] . More involved expressions and fitting strate- 
gies were used in Q. This form of the interaction corre- 
sponds to the class introduced in Eq. (0} and therefore 
we can use Eq. (© to obtain v nn > = J2a (P a )* n (P a ) n ' 
where (P a ) n = j dr<j> n (f)P a (r). For each of the P a (r) 
factors, the corresponding term is nothing but a rank 
one Kronecker approximation to the exact v nn i matrix 
of the interaction. The vectors {P a ) n considered in Refs. 
0~H ITT1 - IT^ are inspired by physical requirements and 
the free parameters entering their definitions can be used 
(as it is the case) to find an optimal (in physical terms) 
approximation to v nn > by the requirements of reproduc- 
ing as well as possible some nuclear properties. From 
our previous numerical study we can conclude that our 
"blind" approximation is also well suited to deal with the 
problem. 
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D. Ultraviolet divergence 

From the previous results one could conclude that, be- 
cause the low rank separable expansion is a kind of "short 
range" expansion, it should show some sort of ultravio- 
let divergence like the pairing matrix elements of a zero 
range contact interaction. In this section I will show that 
this is not the case. First of all, I will remind the reader 
about the origin of the ultraviolet divergence. It is a 
direct consequence of the fact that the nuclear matter 
pairing matrix elements 

- G{k x MWi - n)\k[,k' 2 ) = - —J^ 6(K - K') 

(17) 

of a contact interaction do not depend on the momen- 
tum transfer vector q = k — k! (we have introduced the 
vectors K = -^75(^1 + k 2 ) and k — (kx — k 2 )) favoring 
infinitely high momentum transfers unless some cut off is 
introduced. In the one-dimensional HO case, the struc- 
ture of the general matrix elements is more difficult to 
visualize, but for the sake of discussion, we can restrict 
to the specific matrix elements —G(n,0\8(xi — cc2)|n,0) 
which, as can be shown easily using the preceding for- 
mulas, are independent of n and given by -G/V2. The 
constant value favors, as in the nuclear matter case, the 
scattering of the n — state into high n states. For an 
one-dimensional Gaussian of range \x the expression of 
the previous matrix element is not as simple as in the 
case of the zero range force, and therefore it is preferable 
to plot it as a function of n in order to discuss its prop- 
erties. In Fig. [5] we have plotted the matrix elements 
(n, 0| exp(— (xi — X2) 2 1 '[i 2 )\n, 0) as a function of n as well 
as three of the low rank approximations corresponding 
to Nc = 60 and Lc = 1 and 3 and Nc = 10 and Lc = 3 
(see caption for details). The conclusions from this plot 
are that the matrix element is a decreasing function of n 
quenching the promotion to high n states and therefore 
the ultraviolet divergence. Also, the similarities between 
the exact matrix elements and the two typical low rank 
approximations with Nc — 60 for all values of n are re- 
markable. In the case Nc = 10 and Lc = 3 the approxi- 
mation is reasonable up to n = 15 and from there on the 
matrix elements are too small as compared to the exact 
ones. In any case, we conclude that the low rank ap- 
proximations maintain the characteristic decreasing with 
n of the matrix elements that is needed to prevent the 
ultraviolet divergence. 

E. General interactions 

As discussed previously the case of a general interac- 
tion involves the evaluation of the matrix v nn i and its 
diagonalization. The evaluation of v nn i can be carried 
out numerically using, for instance, Gauss-Hermite inte- 
gration thanks to the presence of the HO wave functions 
in the matrix elements. The resulting formula is also 




5 10 15 20 25 30 
n 

Figure 5: Matrix element i> n ,o,n,o for the one-dimensional 
Gaussian interaction of range [i — 1.4 fm and oscillator length 
parameter 6 =2.0 fm plotted as a function of n. Full curve 
corresponds to the exact result, dashed ones to the approx- 
imation with Nc = 60 and Lc = 1 and Lc = 3 (the one 
closer to the exact result). The dotted curve corresponds to 
the Nc = 10 and Lc = 3 approximation. 

in the form of Eq. (|6]) as can be inferred by using an 
one-dimensional example 

/oo 
dz (p^(z;b)v(V2z)<f n '(z; b) 
-OO 

N H 

i=-N H 
N H 

= D* m (^v{^hzi))D nll (18) 

i=-N H 

where zi and Wi are the nodes and weights of the Gauss- 
Hermite integration, (p n '{ z ) is the reduced HO wave func- 
tion defined above and the coefficients D n i = (p n {zi). 
Also in this case, if v(r) is short range, the quantity 
Wiv(V2biZi) will decrease as 0, increases. In the case of 
a Gaussian interaction this procedure will reproduce the 
exact answer and will provide an alternative (but equiva- 
lent) description of the previous results for the Gaussian. 
Another alternative to this procedure, that lends more 
analytical results, is to represent the interaction to be 
treated by means of its Gauss transform. The Gauss 
transform is a derivative of the more popular Laplace 
transform, where a change of variables allows to express 
given functions as linear combinations of Gaussian with 
different widths. The best known example of this treat- 
ment is that of the Yukawa potential [23fl 

e _Z= 2 r e^^l^dt (19) 
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that can straightforwardly be used to deal with the 
Coulomb interaction by taking the fj, — > limit (see Ref 
[HI for an early application in nuclear physics). The 
advantages of this method are first that the separable 
expansion for the Gaussian is known analytically and 
can be used straight ahead and second that the Gaus- 
sian interaction is separable along the three spatial di- 
rections allowing to treat the problem as three uncorre- 
cted one-dimensional problems one for each spatial di- 
rection. The latter advantage is specially helpful to deal 
with the Coulomb interaction. Given the long range of 
the Coulomb force the number of terms required for an 
accurate separable expansion of a general matrix element 
is expected to be larger than for a gaussian. However, the 
relevant Coulomb matrix elements for the nuclear case 
involve single particle wave functions which are located 
inside the nucleus and therefore explore the Coulomb po- 
tential in the interior of the nucleus. For those matrix el- 
ements, the Coulomb force can be considered as a short 
range interaction with a range of the order of the size of 
the nucleus and therefore the number of separable terms 
required for an accurate representation are expected to 
be much smaller than the ones required for a general ma- 
trix element of the Coulomb interaction. 



study of two specific examples: the contact delta force 
and the one-dimensional Gaussian interaction show that 
the number of terms of the separable expansion to be con- 
sidered can be substantially reduced without affecting too 
much the accuracy of the approximate formula. The sep- 
arable expansion turns out to be a kind of "short range" 
expansion. The issue of ultraviolet divergences inherent 
to any short range expansion is analyzed. The proposed 
separable expansion opens up a new avenue to deal with 
finite range interactions both in the pairing and exchange 
channel. In the former case, the cumbersome definition 
of the "cut off window" is avoided by the implicit finite 
range. The possibility to extend these considerations to 
other finite range interactions like the Coulomb potential 
is also discussed. 
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F. Density dependent forces 

In the applications of the HFB method it is common 
to find terms in the interaction/energy functional that 
are referred to as "Density Dependent" (DD) terms and 
are given by the general expression 

v(fx,f 2 ) = f(r)G(R) (20) 

Usually, f(r) = S(^/2f) and G — p a , that is, the density 
raised to some (usually non integer) power a. Applying 
the ideas developed in this article it is easy to obtain in 
this case 

JV,n 2 ,n 3 ,ru = X] XI ^n^n 2 f L9 kM^^. (21) 
K L 

where we have introduced the coefficients 

M«n 2 = E M r Zl 2 D nL E NK (22) 

Nn 

given in terms of the spectral decomposition of the ma- 
trices / and G with eigenvalues /l and gx , respectively. 
In the common situation where f(f) — 5(\ // 2f) the /l are 
constant and equal to 1/2 3 / 2 and D n ^ = 4> n (0) indepen- 
dent of L. 



III. CONCLUSIONS 

By means of the transformation properties of the HO 
basis to the center of mass and relative coordinate a sep- 
arable expansion of the pairing and exchange matrix ele- 
ments of a general two body interaction is obtained. The 



Appendix A: Talmi-Brody-Moshinsky coefficients in 
one dimension 

The Talmi-Brody-Moshinsky coefficients in one dimen- 
sion are defined as 

Nn 

where X = + x 2 ) and x — -j^{ x i ~ x 2 ) are the 

center of mass and relative coordinate respectively and 
ip n (x) = (^2 n n\b)-^ 2 H n (x/b)exp(-^) ig thg Qne _ 
dimensional HO wave function written in terms of the 
Hermite polynomials. Explicit expression of have 
already being obtained in the literature [20j but we will 
introduce here another derivation which is simple and 
straightforward. We will take advantage of the generat- 
ing function of the one-dimensional HO wave functions 

G(x/b,t) =exp(-i^ + 2^-i 2 ) =J2xn(t)v> n (x) 

n 

with X n(t) = {^b2 n n\) 1 l 2 t n /n\. Given the 
Gaussian form of the generating function we have 
G(x 1 /b,t 1 )G(x 2 /b,t 2 ) = G(X/b,T)G{x/b,t) where T = 
-^{ti +t 2 ) an d t — -^(ti —t 2 )- Expanding in right hand 
side of the last identity Xn(T) and x„(i) in powers of 
tiand £2 and comparing with the same powers in the left 
hand side we finally obtain 

M Nn _ c ( "l!"2! 

m n x n 2 ~ 0„ 1+ „ 2 ,JV+n I 2 ni+T i2 _/V'! n ! I 
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x E(-W N + )( n ) 

where the selection rule is made explicit. 



Using the generating function we can write G(x/b,t) = 
G{x/b',t') exp(-(l - r/ 2 )t 2 ) where 77 = b'/b and t' = rft. 
Equating equal powers in t in the previous expression and 
using Eq. (|B2[) we obtain 



Appendix B: Transformation coefficients 

In this appendix I supply the expression for the 
expansion coefficients of the one-dimensional reduced 
HO wave function (p n {x/b) = exp(^jp)(p n (x/b) = 
(^/n2 n n\b)~' 1, ' 2 H n (x/b) in terms of the same object but 
for a different length scale b' 

- /b'\ 

<p n (x/b) = D nn' i^jMx/tf). (Bl) 

n'=0 ^ ' 

There are many ways to obtain the D coefficients al- 
though probably the more economical one is to use the 
generating function of the Hermite polynomials H n {x/b) 
which is given by (see previous appendix) 

G(x/b,t) = exp f2^i - f 2 ") = ^y^-HJx/b) (B2) 
V b ) n\ 



H n (x/b) = ^(-l)^ 71 '^ 1 ^ 2 A n , n ,H n ,(x/b') 

n' V 2 / ' 



which is an identity that finally yields to 



D nn ,(rj) = A Bl „/(-l)->- -j- ' „ J, - 1 . 

\n'lj 2— (2^)! 

(B3) 

In this expression, the function A„ ; „/ = |(1 + (—1)™+™ ) 
which is one if n and n' have the same parity (even or 
odd) and zero otherwise has been introduced. Notice also 
that D nn i{rf) = if n' > n in agreement with the limits 
of the sum in Eq. (|B1I) . 
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